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ABSTRACT 


This report summarizes the work performed under NASA Grant NGR-40-002-093 
on the theoretical and experimental aspects of single crystal, silicon photo- 
voltaic devices made by forming a grating pattern of p/n junctions on the light 
receiving surface of the base crystal. Based on the general semi-conductor 
equations, a mathematical description is presented for the photovoltaic proper- 
ties of such grating- like structures in a two-dimensional form. The resulting 
second order elliptical equation is solved by computer modeling to give solutions 
for various , reasonable , initial values of bulk resistivity , excess carrier 
concentration, and surface recombination velocity. The validity of the com- 
puter model is established by comparison with p/n devices produced by alloying 
an aluminum grating pattern into the surface of n-type silicon wafers. Current 
voltage characteristics and spectral response curves are presented for cells of 
this type constructed on wafers of different resistivities and orientations. 
Grating type cells are found to have increased spectral response to the high 
energy (blue) photo region of the solar spectrum, and short circuit currents 
equal to or greater than those of commercial cells. Cells produced to date have 
displayed ' open circuit voltages approximately 10 per cent less than those of 


commercial cells. 
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1. INTRODUCTION 

Since the first publication by researchers at Bell Telephone Labora- 
tories describing the invention of the modern silicon Solar Battery in 
1954^, the silicon solar cells have continued to dominate the field of 
photovoltaic energy conversion. During the past 20 years, two important 
developments in cell structure have led to today's high efficiency silicon 
solar cells . 

(2 3) 

The first important development is the shallow junction cell ’ . 

The concept that the collection efficiency is- maximized if a very shallow 

junction is used, taking into account the surface recombination velocity, 

(4) 

was first pointed out by Loferski in 1956 . Earlier in 1954, 

( 56 )' 

Cummerow'' ’ studied the dependence of collection efficiency on junction 
depth, minority carrier diffusion length, optical constant, under the 
assumption that both front and back surfaces have infinite surface recom- 
bination velocities. During the period 1960-1961, several experimental and 

theoretical studies on the variation of the spectral response of cells with 

(7 8 9) 

junction depth were published 5 * . These studies showed that moving the 

junction closer to the surface improved the spectral response of the cell 
in the blue region and hence gave higher efficiencies, particularly in AMO 
(air mass zero) sunlight. 

The second important development is the gridded cell. It was pointed 
out in an early report of Chapin , Fuller and Pearson ^ ^ that series resist- 
ance was a major cause of power loss in the initial cells made at Bell 
Telephone Laboratories. In early 1959, it was recognized that the incor- 
poration of contact grid lines could greatly reduce the series resistance 
caused by the thin diffused layer; Wolf described the benefits derived 
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from contact grids in a paper published in 1960^°^. Grid contacts were 
quickly incorporated into the design of commercial cells 

Coupled with improvement in the base material properties (e.g., long 
minority carrier life time, in the 8 to 10ys range), the combination of 
the grid contacts and shallow junction- depth resulted in an increase of cell 
efficiency from about 6% to about 10% in the period 1959-1961 ^ \ 

There are two inherent disadvantages associated with the shallow junc- 
tion cells that limit the cell performance. 

a) High surface recombination velocity. (Or according to a different 

( 12 ) 

theory developed by Lindmayer and Allison , the recombination. of 
excess charge carriers generated by the blue-violet light is not con- 
trolled by front surface recombination, but by the presence of a "dead 
layer" just beneath the front surface with an extremely short life time. 
For a n on p cell having junction depth of 4,000A, the n layer formed 
by diffusion of phosphorous into the silicon has a heavily damaged layer 
or "dead layer" of about 1,500A.) 

b) Short minority carrier diffusion length in the diffused skin region. • 

Because of the high temperature diffusion process (over 1,000°C) and 

20 -3 

high impurity concentration (varies from about 10 cm at the surface 

15 -3 

to a value of about 10 cm at the junction less than 0.5 pm below the 
surface) incorporated into the skin region, the diffused skin region is char- 
acterized by a high value of surface recombination velocity (or according to ' 

(12 ) o 

Lindmayer and Allison , by the presence of a dead layer of about 1,500A), 

• 5 

and short bulk diffusion length; over 10 cm/sec. and 1 ym, respectively-, as 
compared to values in excess of 'lOOym for bulk diffusion length- in the base 
material and attainable surface recombination velocities of less than 
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3 f 13 14 ) 

10 cm/ sec. for silicon' * (the surface recombination velocity in a 

15 3 

thermally oxidized silicon (4.5 x 10 phosphorous atoms/cm ) structure 
(15) 

studied by Grove was found to be of the order of 5 cm/sec . ) . 

It was these problems of the standard diffused silicon cell that moti- 
vated this investigation of "grating" type cells (see section 2.3). Unlike 
a standard diffused cell which has its front surface covered completely by 
a thin diffused skin region, a grating type cell would have part of its 
front surface covered by a grating of finely spaced fine grid lines which 
form a potential barrier with the base substrate. The grating structure 
junction could be a diffused p-n junction, an alloyed p-n junction or a 
Schottky barrier jucntion. In such a cell, most of the surface exposed to 
the incident light would be the "virgin" silicon constituting the base 
material. It was anticipated that low surface recombination velocities 
would be attainable on this material in addition to the long diffusion 
length in this region (over lOOjim as compared with about lym in the diffused 
skin region), and that the "blue response" of the cell would be increased. 
Its red response would be relatively unaffected provided that the spacing 
between grid lines were smaller than the diffusion length in the base 


material. 


The development of the theory of grating cell, of a fabrication proce- 
dure and the subsequent measurement of the characteristics of the cells 
are the subjects of the work reported here. 
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2. THEORY 


2 .1 General Governing Equations 

In a semiconductor, the behavior of charge carriers, is governed by 

a) Poisson’s Equation, b) Space Charge Equation, c) Drift-Diffusion 
Equation, and d) Continuity Equation, 
a) Poissons' Equation, 


where E is the electric field, 
constant . 



p the charge density. 


and e the dielectric 


b) Space Charge Equation 

p(r) = q(p(r) - n(r) + N + (r) - N_(r)) 
where r is the position vector, p hole density, n free electron density, 
N + (N ) localized, immobile positive (negative) charge density. 

c) Drift-Diffusion Equation, 


<? = qp Ep - qD Vp 

p H p * 


J = qp En + qD ^n 
n u n . H n 

where ^(3^) is the hole (electron) current density, Pp(p R ) the hole 

electron) mobility, D (D ) the hole (electron) diffusion constant. 

p n 

d) Continuity. Equation, 


|E = _ i V -3 - R + G 
at q P P P 

|E= i. v • 3 - R + G 


at 


n n 


n 


where R (R ) is the net rate of recombination for holes (electrons), G (G ) 
p n P n 

is the rate of generation of holes ( electrons ) due to. non -thermal source , 
such as light irradiation, high energy particle bombardment. 

For the steady-state, field free case in nondegenerate extrinsic 
material and under the assumption of low-level injection (i.e., n’, p' « n Q 
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in n-type material , and n * , p 1 « p Q in p-type material , where p * (n ’ ) is 
the excess hole (electron) density, p Q (n Q ) is the thermal equilibrium hole 
(electron) density), R^R^) ma y expressed as p'/x^ (n’/x^)^®^. The 
Drift -Diffusion Equation and the Continuity Equation can be combined into one 
so-called semiconductor diffusion equation, which has the form 


D V 2 p' - 2- + G =0 
P t P 

. in n-type material 

(2.1a)’ 

P 

i 


2 n f 

D V n T - — + G =0 
n x n 

n 

in p-type material 

(2.1b) 


where x (x ) is the minority charge carrier life-time in a n-(p-) type 

hr ** 

material . 

2.2 One-Dimensional (Conventional) p-n Junction Photovoltaic Cell 

In the case of the standard photovoltaic cell, the spectral response 
can be obtained by solving the differential equations in a one-dimensional 
model as shown in Fig. 1. Assume there is no reflection of the incident 
light from the front surface, the absorption of the incident photons with 
energy larger than the energy gap of cell material is goverened by 
Labmert's Law (i.e., equal layer of equal thickness absorbs an equal frac- 
tion of the light which traverses it), and the quantum yield is unity . 

Then the number of electron-hole pairs generated by light with wavelength 
A per unit volume G^tG^) in Eq. 2.1 is aN o (A)e -aX , where a is the optical 
absorption constant of the material and is a function of wavelength, N q (X) 
is the incident light flux density with wavelength A in photons/unit 
area-sec . ) . 

The junction is at a distnace H below the front surface. The thickness 
of the wafer is b, and usually b » £. 
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The one-dimensional form of Eq. 2.1 is 

(2.2a) 
(2.2b) 

Boundary conditions must be specified for the junction and the surfaces 
of the cell. The case of the short-circuited junction (i.e., when there is 
no potential difference appearing across the junction, and the current out- 
put is equal to the light generated current) is considered. For the p/n 
cell (p-type surface, n-type base) shown in Fig. 1, the boundary condition 
at the junction for a short circuit junction has the form, 

n'U) = p’U) = 0 

because according to the assumption of quasi-equilibrium, 
n’(£) = np 0 (e^^^-l) , p’(£) = p no (e^^-l)^^ , and V = 0, where V 
potential difference appearing across the junction, p nQ (n ) is the 
equilibrium minority carrier concentation in n-(p-) material; and at 
surface , 

J n (0) = qs o n'<0) = ,D n f£-| x = Q 

J p (b) = qS b p'(b) = -qD p fj£j x = b 

where S , S, are surface recombination velocities on the front and back 
o’ b 

surfaces respectively. 

Only the case of a thin space charge region is considered here,’ there- 
fore, we can neglect recombination, generation or scattering of the carriers 

in the p-n transition region, in which case Shockley’s diode theory of the 

(18 ) 

p-n junction is valid . The electron and hole contributions to the 


(2.3) 

is the 

thermal 

the 

(2.4a) 

s ■ 

(2.4b) 


D \£i.El +aNe 

.p 3 2 T o 

r X p 

n 3^n' n* . 

D - „ + aN e 

n 3 2 t o 

x n 


-ax 


-ax 


= 0 on the n-side 


= 0 on the p-side 


t 
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total current across the junction are given by: 


J = 
n 


qD 


dn 1 j 
n dx j 


A 


J 

P 



dp 1 | 

dx |a * 


(2.5a) 


(2.5b) 


The ratio of the short-circuit current to the incident photon flux of 
wavelength A neglecting reflection is referred to as the collection 
efficiency Q at wavelength A, 


dr, + d J,, J 

n P _ n p 

^ ~ qN ” qN + qN ^skin * ^base 
o . o o 


( 2 . 6 ) 


Its value specifies the number of electrons flowing through the cell 
per incident photon. 

Solution of the above equations with the appropriate boundary conditions 

yields expressions for the skin and base contributions to the collection 

. (7) 

efficiency : 


. v A -v A 

Q skin * — 2727 <B nV " - W 

a(l-v /a ) 
n 


-aJU 
ae ) 


(2.7) 


where 


-V £ 

3 n = (e n (h + a) - e“ a£ (h + V jp^ A n 


v £ 


y = (e a£ (h - v )-e n (h + a)*)/A 
n n n 


v A -v A 

= e n (h + v r ) - e n (h - v n ) 


h = S/D 
o n 


v = (D T )“ 1/2 = L 1 , 
n n n n ’ 


and 


-ba 


(-6 v e P 


v (b-A) 


-v (b-A) 


base 2. 2, P pp 

o(l-v p /a ) v v 


+ y v e 
P P 
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where 

-v(b-£) ,, , 

6 p = (e p (h-a) - e a(b ~ £) (h+v p ))/A p 

,, „s v (b-£) 

Y p = (e“ A \h-v p ) - e P (h-a))/A p 

v (b-A) -v_(b-£) 

A = e P (h+v ) - e p (h-y ) . 

P P P 

h = V D p 


v = (D t ) 1/2 = L 1 . 
P P P P 


To obtain the I (short-circuit current) when a cell is illuminated 

SO 

by a broad band source, this calculation must be repeated to obtain a 
collection efficiency value for each wavelength used for the illumination • 
(because of variations in optical obsorption constant with wavelength). 

An integration must then be performed to account for the dependence of 
illumination intensity on wavelength. 

2.3 Two-Dimensional Grating-Type Diffused p-n Junction Photovoltaic Cell 

The diffused grating type cells are cells whose light receiving surface 
is not covered completely by a diffused junction; the surface is partially 
covered by a diffused grating or grids of collecting junctions. The top 
view of a grating-type cell is shown in Fig. 2. 

The diffused grating lines have a width "2a" and the spacing between 
the centers of these lines is "2b". Because of the periodicity of the 
grating structure, the collection efficiency of a grating cell can be ob- 
tained by solving the semiconductor diffusion equations in a unit of width 
"b" as shown in Fig. 3 and also because there is little variation along the 
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direction of grid lines a two-dimensional model is adopted. (Note that 
here we use n-type material as base, because n-type material is more suit- 
able for making alloyed p-n junction and Schottky barrier cells). The 
governing equations and the conditions to be satisfied on the boundaries of 
the unit cell are shown in this figure. Because the neighboring cells are 
identical, due to symmetry requirement, we have to have 9n ' /9y , 3p ' /3y 
equal to zero at the boundaries between neighboring unit cells. 

The two-dimensional diffusion equation. 


o -ax 
r— e . = 0 


where U can be either n' or p’, does not have a closed form solution as in 

the one-dimensional case. It is a second order elliptical equation which 

(19 20) 

can be solved approximately by finite-difference techniques * which 

replace the continuum by a discrete model, thus reducing the problem to a 

system of simultaneous linear difference equations. 

The matrix of the resulting matrix equation is real, symmetric, with 

positive diagonal entries and non-positive off-diagonal entries, and is 

irreducib'ly diagonally dominant (the same thing is .true for one-dimensional 

(19) 

diffusion equation). In addition to these properties, it is also a 

sparse matrix, i.e.., a large percentage of the entries of this matrix are , 
zero , therefore , the problem can be more appropriately solved by iterative 
methods. Of all the -iterative schemes for solving linear systems, the 
successive over-relaxation method has been one of the most -widely used. 

Its simplicity and firm mathematical basis together with numerical experi- 
ence accumulated over the years make it a convenient, efficient, and reli- 
(19 20) 

able method. * Therefore, the successive over -relaxation iterative 
method was chosen for the computation. (See Appendix) 
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In order to test the effectiveness of the computational method (suc- 
cessive over -relaxation iterative method)-, the response- for a cell completely 
covered by a diffused region (one-dimensional case) was computed and the 
results were compared with -the closed form solution which exists for 
that . case . 

Figure 4 compares the base responses based on the exact solution and 

on a numerical approximate solution for the one -dimensional case. The 

case S 2 = «° corresponds to the case of a metal contact covering the whole 

back surface. The curves show that the numerical approximate solution is 

3 —1 

very good for a larger than 2 x 10 cm , however , the ' numerical approxima- 
tion underestmates the low a response. The profiles of the excess carriers 

concentration per incident photon with a as a parameter, in the base 

* 

region from the exact solution are shown in Fig. 5. We can see that as 
a keeps increasing, the slope at the two ends of the profile increases first 
and then decreases, but we used the same discretilization scheme for all a 
in the numerical computation, that means this approximation scheme is more 
accurate for a with smaller slopes, and indeed. Fig. 4 reflects this. 


* 


P’(x) 


aN e 
o 


-aZ 



b-& 

-a(b-Jt) L p 
e - e r 


b-i b-Jt 



e 


L P e Lp - e~ a(b ~ S ' ) 


b-Z b-Z 



X 

L 

P 


+ e 


-ax 


• e 


for & $ x * b 
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Figure 6 compares the skin responses based on the exact solution and 
on a numerical approximation for the one-dimensional case. In this case 
the numerical solution deviates from the exact solution in the high o region. 

The reason for this is simply that as a increases, the variations of the 
profile of excess carriers in the diffused skin region gets, larger and with 
the same discretilization scheme, the approximation becomes less accurate. 

If we use finer mesh spacing, the deviation could be smaller. Note that, 
in a grating cell, the fraction of the cell area covered by the grating 
collectors is a/b; consequently, thd contribution to the spectral response 
from the diffused skin region is about a/b times of that for one-dimen- 
sional case and the error in total response is accordingly reduced. 

These comparisons of the exact and numerical solution for the one-di- 
mensional case provided an estimate of the reliability of calculated spec- 
tral response curves of grating cells. In the computation for diffused 
grating cells, the sample thickness was set at 100 pm; the junction depth 
was set at 1 pm; the minority diffusion length in the base — = 300pm, 
in the skin = 1 pm; the diffusion constant in the base and in the skin 

O 

D.. = D , . = 10 cm /sec. (for a p/n n-type cell, this corresponds a doping 

base skin 

15 _ 3 

density 10 ' L cm in base material and 10 cm in skin material). The 

surface recombination velocity of the diffused region was assigned value 
5 

= 10 cm/sec , , the surface recombination velocity of the base regxon was 

q 

assigned values S 2 = 0 and S 2 = 10 cm/sec. The number of points imposed for 
discretilization was limited by the cost of the numerical computation, but 

I 

we believe the computation results have supplied us with the maximum information 
needed. The main purpose of the computation was to determine how the overall 
response of a grating cell changes when the a/b ratio changes and how the 
response compares to that of a conventional diffused solar cell. 
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Figure 7 shows the results of the computation for b/a ratio = 5, a = 2pm, 

and the indicated values of surface recombination velocity Sg. The figure 

shows base response, the diffused grating (skin) response, and total response. 

3 

It is reasonable to expect that if the value of is between 0 and 10 cm/sec, 
the corresponding curves should also fall between the curves for S 2 = 0 and 

3 

S 2 = 10 cm/ sec. 

Figure 8 shows the response curves for cases in which b/a is kept at 
5, but "a l! is assigned two different values, viz., a = 2 and 10pm. The 
diffused grating (skin) responses for the two cases are essentially the same. 
The base response for a = 10pm is overall lower than that for a = 2pm for all a 
Figure 9 shows the qualitative description of the response behavior for 
fixed b/a ratio and different values of a, when all other parameters which 
affect the response are held constant. It shows that as "a M increases, the 
skin response remains the same, but base response decreases monotonically , 
for all o. Note that in this figure, the curves labeled 1 and 2 are taken 
from Figure 8, curve 3 is extropolated to lie between these two computed 
curves. Since b/a is fixed, the ratio between the number of hole-electron 
pairs generated by photons of energy e (corresponds to an optical absorption 
constant ct(e) of the material) in the base region and those in the diffused 
grating (skin) region is also kept constant. However, the average distance 
carriers travel to the junction is essentially the same for the light generated 
minority carriers in the skin region and increases for those in the base region 
In other words, the probability of recombination before reaching the junction 
remains the same for light generated minority carriers in the diffused skin 
region and increases for those in the base region, as "a ,! increases. Therefore 

the skin response should remain the same, but base response should decrease 

* 

monotonically for all a as "a ,! increases. 
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Figure 10 shows the qualitative description of the response behavior 
for fixed "a" and different "b/a" ratios. The skin response is roughly 
reduced to a value of "a/b" times the response for the corresponding one- 
dimensional case (b=a), because the number of minority carriers generated in 
the diffused grating (skin) region per incident photon, per units area (total 
exposed area) is reduced to a/b of the one-dimensionai case, and the average 
distance for these carriers to travel to the junction remains the same. 

In the base region, as ’’b" becomes larger, the average distance 
for minority carriers to travel to the junction (generated by monochromatic 
light) increases, and this is true for carriers generated by monochromatic 
light of all wavelength in the useful range. However, there are two cases 
to consider, 

a) for high energy photons, those for which a is high in the material, 
the number of minority carriers generated per unit area (total exposed 
area) increases rapidly (because when b=a, almost all carriers generated are 
in the skin, very few in the base, and when b=2a, suddenly more than half 
of these carriers generated are in the base); 

b) for low energy photons, those for which a is low in the material, 
the number of minority carriers generated per unit area remains the same 
(because even when b=a, almost all of these carriers generated are in the base 
region). Therefore, the base response increases for high energy photons and 
decreases for low energy photons as b/a begins to increase from unity. 

Eventually, however, as "b" continues to increase, the average number of minority 

i 

carriers generated per unit area by high energy photons approaches a constant 
value, and no longer is enough to compensate for the effect (more recombinations) 
caused by the increase of the average distance minority carriers must travel 
to a collecting junction. Thereafter, the base response Starts to drop for 
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all ct, as indicated in Fig. 10 for the case b>>a. 

Note that in Fig. 10, the curve with b=b 2 is reproduced from Fig. 7 
for the case S 2 = 0, the curve with b = b Q is the exact solution for the 
corresponding one -dimensional case. The other curves have been plotted 
qualitatively according to the reasoning above. 

In order to compare the overall collection efficiency of a diffused 
grating and a corresponding one -dimensional cells under AMO illumination , 
the short circuit current for both cells are calculated from the integral : 

I op = q I 00 n(e)N (e)de (2.10) 

SC Jl.l2eV e 

where 1.12eV is the energy gap of single crystal silicon at room temperature, 

N £ (e) is the photon flux density in AMO spectrum e is the photon energy, 

n(e) is the cell response to a photon of energy e. N £ (e) as shown in Fig. 11 

.( 21 ) 

was computed by Wysocki from Johnson's data for the spectral irradiance 
at AMO. 

For convenience in calculation, the Johnson curve of Fig. 11 was trans- 
formed into a curve of incident flux density as a function of wavelength X, 
through the following relation : 

N. (X)dX = -N (e)de (2.11a) 

A e 

V x) . -n e <o ft - ( E ) f- (2.iib) 

A 

the resultant curve is shown in Fig. 12. 

The response vs. photon wavelength for both the grating cell and the 
corresponding one -dimensional cell is shown in Fig. 13 which was transformed 
from Fig. 14 with the help of a vs, hv curve (Fig. 15)^ 22,23 \ Figure 14 
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shows the total response vs. a curves of the diffused grating cell whose 
response was calculated and displayed in Fig. 7 (the case b/a = 5, a = 2pm, 

= 0, S 2 is the surface recombination velocity of the- "virgin" silicon 
between the grating lines, and a zero recombination velocity on the back 
surface), and the corresponding one-dimensional cell. By combining the 
information exhibited in Fig. 12 and Fig. 13, the short circuit current, from 
these two cells exposed to AMO sunlight can be computed by using the relation. 




1.11pm 

0 


n(X)N x U) dX 


( 2 . 12 ) 


The ratio of short circuits between the diffused grating cell and the 

corresponding one-dimensional cell is roughly 1.04. The response curve used 

for the diffused grating cell may seem to be a very favorable case , since 

in computing this curve, it is assumed that the surface recombination velocity 

= 0 for the "virgin" area, i.e., the area between grating lines. However, 

the grating structure used in the computation is not an optimized one and Fig. 7 

3 

shows that there is not much difference between curves with S 2 = 10 cm/sec. 
and S 2 = 0, and the value of S 2 attainable is expected to be much less than 

q 

10 cm/sec. (The surface recombination velocity in a thermally oxidized silicon 

(15 ) 

structure studied by Grove was found to be of the order of 5 cm/sec.). 

The assumption that the back surfaces of these two cells have zero 
recombination velocity is for the simplication of computation, however, the 
order of magnitudes of these two short circuit currents is not expected to 
change if an infinite recombination velocity on back surfaces is used. 

Because for those carriers generated near the back-surface, they are dominated 
by the condition on the back surface, and as long as the optical absorption 
constant "a" does not change due to the incorporation of diffused grating 
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structure, the front structure of the cell would have little influence on those 
carriers generated near the back surface. 

In summary, these calculations confirm that it is possible to design 
grating cells with higher short-circuit current circuit under AMO illumination, 
whose response at long wavelength is less than that of a corresponding 
conventional cell, but whose response at blue-ultraviolet region is larger. 

For a given conventional (one dimensional) cell, and a given minimum admissible 
value of "a", there are an infinite number of corresponding grating cells 
(in the sense that when "b" = "a", these grating cells become identical to 
the one-dimensional cell), a cell among all these cells with the highest I SC 
under some certain illumination spectrum is called the optimized grating cell 
voider that illumination. The optimized value of "b" would be strongly affected 
by the minority carrier diffusion length in the base region; for an exceedingly 
short base diffusion length, the optimized "b" may be even equal to "a". 

2.4 Two-dimensional grating-type p-n junction photovoltaic cell with opaque 

grating lines 

The structure of a grating type p-n junction photovoltaic cell with opaque 
grating lines is like that for a diffused grating type p-n junction cell, but 
with an opaque metal grating overlapping the grating structure p-n junction 

to form an ohmic contact; the junction can be either alloyed type or diffused 
type. 

Assume that the metal grating is completely opaque, then the response 
comes only from the base region. For the case, "b"="a" (the cell is completely 
covered by the metal), the response vanishes. For the case "b">>"a 1 ’ and 
"a"<<L^ ase , the grating intercepts only a small fraction of the light incident 
on the cell and the response is essentially the same as that of a diffused 
grating cell having the same grating structure. For intermediate values of 
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the b/a ratio and a « , the high a response of the base region is the 

same for the two cases , since for both cases it arises from absorption of 
photons in the near surface region of the base silicon between grating 
lines ; the low a base response is approximately (b-a )/b times the diffused 
grating cell value, since it is assumed that the light incident normal to 
the surface continues to travel along this direction without deviation from 
its path, and the grating "shades" a/b of the exposed area. Because the light 
may not be normally incident and because of scattering the actual response 
of real opaque line grating cells for low a photons should lie between the 
dashed and the solid curves shown in Fig. 16 and Fig. 17. 

Figure 16 is a qualitative description of the base spectral response 
behavior for opaque line grating cells with a fixed b/a ratio and different 
values of "a", which is obtained by modifying the base spectral responses of 
corresponding diffused p-n junction grating cells as shown in Fig. 9. Note 
that for low energy photons , the base response is approximately (b-a )/b 
of the value of that of a corresponding diffused p-n junction grating cell. 

In this figure, we have assumed b/a = 5. 

Figure 17 is a qualitative description of the base spectral response 
behavior for opaque line grating cells with a fixed value of "a" and different 
values of b/a ratio, which is obtained by modifying the base spectral 
response of corresponding diffused p-n junction grating cells as shown in 
Fig. 10. Note that for low energy photons, the base response is approximately 
(b-a)/b of the value of that of a corresponding diffused p-n junction grating 
cell. In this figure, we have assumed b Q = a, b-^ = 2a, b 2 = 5a, bg » 

From Figure 16, it is concluded that for a fixed "b/a" ratio, the 
response at all a decreases as "a" increases, and referring to Fig. 15 
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"absorption spectra of single crystal silicon at 300°K", it is equivalent 
to say that the response at all wavelength decreases as "a" increases. 

Therefore, under the illumination of some specific light source, if "b/a" 
ratio keeps constant, the short circuit current will decrease monotonically 
as "a" increases. 

From Fig. 17, it is concluded that for fixed "a", the response at all 
a increases (referring to Fig. 15, it is equivalent to say that the response 
at all wavelength increases) as "b/a" begins to increase from unity, it 
will reach a maximum value at some certain value of "b", thereafter, the overall 
response will drop as "b" continues to increase, in contrast to the more 
complicated behavior for a diffused p-n junction grating cell as shown in the 
same figure . Therefore, for opaque line grating cells, if "a’’ is kept constant, 
the short circuit current will increase as "b/a" begins to increase from unity, 
it will reach a maximum value at some certain value of "b" (for a minimum 
admissible width of "a", this value of "b" depends on the minority carrier 
diffusion length, and on the surface recombination velocity of the base material), 
and will drop as "b" continues to increase. 

2.5 Concerning possible methods for fabrication of grating type cells 

A number of possible configurations of grating cells were considered. 

The diffused grating p-n junction cell can be made by using microelectronic 
techniques to define a diffused grating p-n junction pattern on a base wafer. 

The steps involved would be as follows. The surface of a n-(or p-) type Si 
single crystal wafer would be covered by a masking SiO^ layer. Photoresist ing 
techniques would be used to define and etch the desired grating pattern into 
the oxide layer. The sample would be subjected to the standard diffusion 
process and a p-n junction grating pattern would be formed on the surface. 
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The wafer would again be covered by a masking SiO x layer. Windows would be 

opened up in the SiO^ layer for deposition of ohmic constants to the diffused 

grating structure. These ohmic contacts could be fine metal lines, narrower 

than the diffused grating lines but forming a contact to every one of the 

diffused grating lines along its length. This method of fabrication of 

(24 ) 

grating cells was attempted by others in our group ; -we have concentrated 
on a simpler system, namely the alloyed-aluminum grating cell. The p-n 
junction grating pattern is formed not by diffusion, but rather by alloying 
A1 into silicon. The aluminum network also constitutes the ohmic contact 
to the grating. Consequently the number of steps required to produce a 
grating cell is reduced. As we show below, it is even possible to omit 
the oxidation steps required to form a mask for the diffusion process. 

The aluminum alloying process requires heating the wafer up to the vicinity 
of 700°C for a few minutes; it would be preferable if this heating step could 
be eliminated or at least if the temperature were reduced. A Schottley 
barrier type grating cell would make this possible. The reverse saturation 
current I q of a barrier type photovoltaic cell plays a very important role 
in setting the efficiency of a cell. It can be seen from a study of the 
simplified current -voltage characteristics of a barrier-type photovoltaic 
cell. The relationships between load current I and load voltage V has the 
form: 

I = I ( e qV/AkT - 1) - I (2.13) 

O bL 

where I q is the reverse saturation current , Ig c is the light generated 
current: A is a constant whose value would be unity in a classical Schottley 

i < 

p-n junction but whose value in barriers formed on silicon lies between 1 
and 2 . In this expression the effects of internal series resistance R and 

u 
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shunt resistance R ^ are neglected. From this equation, it follows that 
the open circuit voltage V qc is given by 


„ AkT . , sc , , . 

V = — — lnC-^jr- +1) 

oc q 1 

u o 

The load voltage at the maximum power point V is in turn given by 

qV " qV I. 

(1 + ^ ) exp ^ + 1 


(2,14). 


AkT 


(2.15) 


High efficiencies require high values of V which in turn require low value 
•of I q . The value of I gc should not depend on barrier height, provided that 
some sort of potential barrier is present. The value of I Q , however, is 
extremely sensitive to barrier height. In the case of the thermionic emission 
model v 7 of a Schottky barrier junction 


< 1*8 

I o « exp (- ) (2.16) 

Therefore, low values of I Q require high value of barrier height 4> ^ • The 

(25 ) 

value of <j>g for gold on n-Si is about 0.8/eV v , Au is a "lifetime killer" 

in Si, and therefore its employment in Si grating cells is undesirable in 

( 26 ) 

photovoltaic application. The value of <}> D for A1 on n-Si is about 0.76 eV v , 

p 

however, a literature research showed that the reverse saturation current I 

o 

for Au-n-Si Schottky barrier is significantly (several order) higher than 
that achieved in Si p-n junctions 5 ^ \ f For Al-n-Si Schottky barrier 
junction, the I is even higher, it was for this reason that we chose to 
concentrate on alloyed aluminum p-n junction grating cells. 
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3. EXPERIMENT 

3.1 Fabrication of aluminum-n-silicon alloy p-n junction cells 

If an aluminum film is placed in contact with n-type silicon, and the 

resultant metal-semiconductor system is heated to its eutectic temperature 

(576°C) in a reducing atmosphere, there forms at the interface of the two 

materials a thin melt having the eutectic composition 11% silicon and 89% 

(29) 

aluminum . At this temperature , the molten region continues to enlarge , 

with the two components dissolved in the liquid phase at the eutectic composi 

tion, until the aluminum is exhausted. The aluminum is the limiting element 

because of the thinness of the aluminum film (less than 1pm) compared to 

the thickness of the silicon wafer (~250 pm), and the high percentage of 

aluminum in the eutectic composition. 

As the aluminum-silicon system is heated above 576°C, the aluminum- 

silicon phase diagram (Fig. 18) shows that the percentage of silicon in 

the liquid phase increases. When the temperature of the system is lowered 

through the eutectic point, excess silicon in the liquid phase is rejected 

and a regrowth region forms at the liquid-silicon interface. The silicon 

base region acts as a seed for recrystallization of the rejected liquid 

silicon, and the regrown region is a single crystal with the same orientation 

(30 31) 

as the undisturbed base region * . The phase diagram does not show 

it , but there is a small amount of aluminum which remains dissolved in 

18 3 ( 32 ) 

silicon in the solid phase, approximately 7x10 atoms/cm . Since 

the donor impurity concentrations in resistivities commonly used for 

solar cell manufacture, i.e., 100 ohm-cm to 0.01 ohm-cm, are in the range 
14 18 3 

10 to 5x10 atoms/cm , the regrowth region will be converted into p-type 
and an abrupt p-n junction will be formed at the regrowth region -base wafer 


interface . 
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The following steps have been followed in making Al-n-Si alloy p-n 
junction grating cells, 

a) The prepolished n-silicon wafer having resistivity in the range 

0.1 to 10 ohm-cm is cleaned in accordance with standard procedures commonly 
used for fabrication of microelectronic circuits (i) The wafer is cleaned 
ultrasonically in trichloroethylene for 30 sec . , followed by cleaning in 
acetone for 30 sec. and in methanol for 30 sec. (ii) It is washed in a 
1:1:5 solution of NH 4 0H, and deionized H 2 0 at 80°C for 15 min. 

(iii) It is washed in a 1:1:6 solution of HCL, ^ 2° 2 an< ^ deionized H 2 0 at 
80°C for 15 min. (iv) It is rinsed in deionized H 2 0. (v) It is dipped in 
buffered HF for 20 sec . ) . 

o 

b) A thin layer of aluminum (about 7,000A) is deposited on the polished 

—6 

surface by vacuum evaporation at a pressure of about 5x10 torr. 

o 

c) Two thin layers, the first of titanium of a few hundred A and the 
second of silver of about 1pm are deposited on the other side of the silicon 
wafer . 

d) A layer of A-Z positive photoresist (AZ-111, Shipley Company, Inc., 
Newton , Mass . ) is spread over the aluminum covered surface . 

e) The photoresist is exposed through a mask which defines the grating 
pattern or patterns. 

f) The photoresist is developed with A-Z developer. It is found that 
the developer is alkaline based which removes both the photoresist and the 
aluminum from the regions which are intended to be the open spaces between 
the lines in the grating pattern. This eliminates the metal etch step. 

g) The remaining photoresist is removed by dissolving in acetone, thus 
leaving behind on the wafer an aluminum pattern of the desired form. 

h) The wafer is heated in a reducing atmosphere at temperature in the 
vicinity of or above, the silicon-aluminum eutectic temperature (576°C). 
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During the heating process the aluminum alloys with the silicon forming 
p-n junctions and the Ti-Ag layer forms an ohmic contact to the base. 

i) Leads are attached by bonding Cu wires with "silver print" paint to 
the A1 on the light receiving surface of the cell. 

In order to explore the relationship between open circuit voltage V q 
and heating process, a group of samples having the configuration shown in 
Fig. 19 were prepared. Twenty-one samples were prepared from each silicon 
wafer having a resistivity of 2-3.5 ohm-cm and a (100) surface. 

After the pattern had been defined (as described above), the samples 
were subjected to two kinds of alloying schedules. In one of these the 
samples were heated "slowly", in the other, "rapidly". In the "rapid" 
heating experiments, the controller of the furnace was set at the 
desired temperature and the furnace was turned on. It was found that with 
the combination of furnace and controller used, the interior of the furnace 
could reach a temperature in the vicinity of 700°C in about five minutes. 

In the "slow" heating experiments, the controller setting was gradually 

i * 

increased up to the desired temperature; about 15 minutes were consumed in 
attaining 700° C. 

In both cases , the steps in the alloying procedure were as follows . 

The sample was placed inside a quartz tube which passed through a tubular 
furnace. The sample was positioned so that it would be in the middle of the 
tubular furnace. The tube was sealed and was purged with' N 2 gas for 15 
minutes followed by H 2 for 15 minutes before the furnace was turned on. 

The H 2 gas continued to flow through the tube during the heating cycle 
until the temperature in the furnace had returned to room temperature . 

In both the rapid and slow heating cycles, the furnace was shut off 
as soon as the furnace had reached the desired processing temperature, i.e., 
the sample was kept at its maximum temperature for a few minutes only. 



-24- 


The open circuit voltage V was measured by exposing the cells to light 

OC i 

from a tungsten iodide source coupled with a water filter intended to reduce 

the infrared content of the incident light. The intensity was about 
2 

100 mW/cm , the samples were measured in succession at the same position. 

The measurement was performed quickly to minimize temperature increases. 

The aluminum covered "squares" of Fig. 19 are opaque, a photovoltage is 
generated because carriers produced around the periphery of the junction 
can be collected. 

It was found that the values of V qc measured on junctions produced by 

the "rapid" heating process were higher than the values measured on junctions 

produced by the "slow" heating process. Figure 20 shows the measured value 

of V qc for junctions produced by the "rapid" heating process. The curve 

exhibits a broad maximum with V values in the vicinity of 500 mV for 

oc 

alloying temperature between 650 and 750°C. 

3.2 Optimization of grating geometry 
a) Cell design parameters 

The response of the grating type cells was discussed in the "theory" 
sections. The analysis given in section 2.4 predicts that, for an opaque 
line grating cell, when the line width "2a" is fixed, the short circuit 
current I will increase as b/a begins to increase from unity, and will 
reach a maximum value for certain line width "2b ” (this value of 2b should 
depend on the value of 2a, and on the diffusion length and surface recombina- 
tion velocity of the base material). It also predicts that when b/a ratio 
is fixed, I of the opaque line grating cell should decrease monotonically 

SC 


as T, a n increases* 
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In order to verify these points, alloyed Al-n-Si grating cells of two 

different configurations were constructed, the dimension of the individual 

cells and their geometrical arrangement on the base wafer is shown in Fig. 21. 

In configuration 1, cells with a fixed line width "2a" and varied spacings 

were fabricated. The line width was set at the minimum line width which 

can be generated at Brown's Microelectronic Lab (about 5ym) and the spacing 

was varied between 30 and 250 pm. The area covered by the grating pattern 
2 

were 50x50 mils • Six groups of three cells each were made on five n-Si 
wafers having different resistivities and orientation, viz. } o.l ~ 0.3 
ohm-cm <111>, 2 ~ 3.5 ohm-cm <100>, 5 ohm-cm <111>, 10 ohm-cm <111> , and 
8 r 12 ohm-cm <100>. These wafers are designated as wafer I, wafer II, 
wafer III, wafer IV and wafer V respectively. The actual values of "2a" 
and "2b" are given in Table 1. Figure 2 is a picture of the top view of 
group 2 cells. 

In configuration 2, cells with a fixed line spacing to line width 

(b/a) ratio and six different line width "2a" were constructed. The b/a 

ratio chosen is 3, and the value of "2a" ranges from 10 to 80 ym. Six groups 

of three cells each and three groups of 2 cells each were made on a n-Si 

wafer with 1 ~ 2 ohm-cm resistivity and <100> or ientation , ’ and this wafer 

.2 

is designated as wafer VI. Group 9 consisted of a two 50x50 mils "picture" 
frame", i.e., no grating pattern covered the surface ; they were used to 
estimate the effect caused by the boundary of small area cell. Group 7 and 
Group 8 have the same grating structure as Group 2 and Group 4, respectively; 
however, the former also have four grating lines at right angle to the main 
grating pattern. The actual value of "2a" the line width and "2b" the line 
spacing of these groups are given in Table 2. 
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All these cells were fabricated according to the "rapid" heating cycle 
described before . 

b) Measurement of I 

sc 

The short cirucit current I of these cells were measured for each of 

sc 

these cells exposed to a tungsten iodide lamp. The light intensity was 

2 

set at a value which resulted in a short-circuit current density of 26.3mA/cm 
2 

from a 2-cm commercial silicon solar cell without anti-reflection coating. 

This value was chosen because the same commerical cell generated a short- • 

2 • . 

circuit current density of 26. 3mA/ cm when it was exposed to an AMO solar 
similator. The highest values of short circuit current for each of these 
groups on wafer I through wafer VI are recorded in Table I and Table. 2. 

In Table 1, it is observed that wafer V and wafer IV have almost equal 
I for corresponding (same "a", same "b") grating cells, this implies 

SO 

that their minority carrier diffusion lengths are almost identical. 

The implication is drawn because of the following reasons: 

i) all the wafers were subjected to the same fabrication process, thus 
they are supposed to have identical surface condition and therefore the 
surface recombination velocity of the base material, 

ii) all the wafers have a thickness of 250 ym ± 50 ym, -in this range 
the cell thickness have little effect on collection efficiency; therefore 
the only major factor left that would affect the collection efficiency 
and Consequently the I of an opaque line grating cell is the minority 

SC 

carrier diffusion length in the base material. 

Hence, by comparing the I for corresponding grating .cells in these 
wafers, it appears that the values of the diffusion length in the base 
material would vary according to the following descending order: wafer V, 
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wafer IV, wafer II, wafer III and wafer I, since this is the descending 
order of magnitude for the I of corresponding grating cells in these 
wafers . 

From Table 1, it is found that for 2a~ 5 urn, -the optimum line spacing 

"2b" for wafer V and wafer IV is i 250 pm, for wafer II, it is about 

127 pm, for wafer III, it is about 64 pm, and for wafer I, £ 32 pm. 

Table 2 lists the I 's for cells fabricated on wafer VI, with a fixed 

sc 

b/a ration =3. It is evident from this table that as "a" increases, 

I decreases monotonically. The only exception to this monotonic decrease 
is group 3, we attribute this deviation from the pattern to unspecified 
imperfections in the cells of this group. We also expect that the rate of 

t 

decrease of I with increasing "a" would have been faster than what showed 

in Table 2. This might be caused by the interference of the side lines of 

these grating cells. The length of each of these four side lines and grating 

2 

lines of these 50x50 mils grating cells is only 1,270 pm, and the line 
spacings have lengths as long as 127 pm and 254 pm; therefore the strictly 
two dimensional model is somewhat violated for these large spacing grating 
cells . 

c) Spectral response measurement 

A Jarrel-Ash monochromater and tungsten-iodide source were used to make 
spectral response measurement for' all cells on the six wafers over the range 

o o 

4,000A ~ 10,000A. The output of each cell tested was measured as a function 
of the wavelength of incident radiation and relative to the output of a 
standard cell whose spectral response was measured at NASA Lewis Research 
Center*. The results were plotted as response per photon and were normalized 

* 

The author wishes to thank Dr. Henry Brandhorst and Mr. Russell Hart of NASA 
Lewis Research Center for arrangement and assistance in using their facilities. 
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to unity. at peak response for each individual cell. All these cells have 
a- common feature, viz., they are blue shifted cells. Figure 22 shows a 
representative case. It shows a normalized spectral response curve for a 
commercial solar cell with its antireflection coating removed and the 
normalized spectral responses for these cells' on wafer II. All these small 
cells normalized responses are confined within two small horizontal bars 
as shown in the figure. It is evident that these grating cells have higher 
blue response and lower red response than that of a commerical silicon 
solar cell. 

Since the normalized spectral response curves for individual cells on 
each wafer in general have the same .shape, the relative spectral' responses 
for cells on each wafer can be plotted by using their normalized spectral 
response curves and their short circuit currents as listed in Table 1 and 
Table 2. Typical results are shown in Figs. 23, 24, 25, and 26. Figure 23 
shows the case of wafer I. For cells on wafer I, the spectral response' 
decreases monotonically with- increasing "b' ! for b 5 32' yin • this means the 
optimal value of "b' ! for wafer I is £ 32 ym. We expect that for 
b i 191 ym, the spectral response would have continued to .decrease -and not 
reaching some limiting value as shown in the figure. The reason is because 
of the very short diffusion length of this base wafer (see discussions 
of section 3.2b), only those carriers generated near grating lines are 
contributing to the response, therefore the contribution to the response 
by side lines is increasingly important for cells with large value of ' b“ 

(for b = 254 ym, there are only 4 grating line collectors in a 50x50 mils cell). 
Figure 24 shows the case of wafer III. For cells on. wafer III, the overall 
response reaches a maximum when "b" is about 64 ym, and- decreases monotonically 
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as "b" continues to increase. Figure 25 shows the case of wafer V. For 
cells on wafer V , the overall response keeps increasing as "b" increases 
from 32 pm, the optimal value of "b" for this wafer is 5 254 um. Figure 26 
shows the case of wafer VI, for which b/a is fixed at about 3. For cells 
on this wafer, the overall response decreases monotonically as "a" increases. 

In summary, these results confirm the qualitatively description shown 
in Fig. 16 and Fig. 17 (via a relation between a and A shown in Fig. 15), 
and the statements in section 2.4, viz., i) for fixed b/a ratio, the overall 
spectral response decreases as "a" increases, ii) for fixed "a", the overall 
response increases as b/a begins to increase from unity, it will reach a 
maximum value at certain value of "b" (for wafer 1, this b is 5 32 pm, for 
wafer III, this b is about 64 pm, for wafer V, this b is i 254 pm), there- 
after, the overall response will drop as "b" continues to increase. 

2 

3.3 Fabrication and measurement of 1-cm grating cells 
2 

Four 1-cm grating cells were fabricated on a 2 ~ 3.5 ohm-cm, (100) 

plane of a n-silicon wafer by using the "repid" heating cycle described 

above. The grating pattern consists of 7.5 pm grating lines spaced 127.2 pm 

apart. Its magnified picture is shown in Fig. 27. 

The outputs of these cells were measured under a simulated AMO Sunlight 

2 

at NASA Lewis Research Center. The intensity was set at 136 mW/cm and the 

spectral distribution is a good approximation to the AMO spectra. The I-V 
2 

curve of a 1-cm grating cell is shown in Fig. 28, which also includes 

2 

the I-V curve of a 2-cm commerical Si cell whose antireflection coating has 

2 

been removed. (The current reading for the 2-cm cell has been adjusted 
2 

to current / cm ) . 
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2 

The following point can be made. The 1-cm grating cell produced a short 

circuit current comparable to and indeed a little greater than that produced . 

by the commercial cell, (with an antireflection coating, both cells should 

produce I about 1.4 times the uncoated cell values). However, the open 

circuit voltage of the grating cell (0.47 V) is smaller than that generated 

by the diffused commercial cell (0.55V). The maximum power point for this 
2 

1-cm grating cell occurs at 0.34 volts and 22.5 mA, this corresponds to a 

measured maximum efficiency of about 5.6%. The commercial cell of Fig. 28 

has a measured maximum efficiency of about 8% with its AR (antireflection) 

coating removed, it occurs at 0.45 volts and 24 mA. The addition of an AR 

coating usually increases output by a factor of 1.4. Thus the commerical 

cell with AR coating would be expected to have an efficiency of about 11.2% 

2 

and the 1-cm grating cell with AR coating would be expected to have an AMO 

efficiency of about 7.9%. As we noted above, the short circuit currents 

of these two cells are comparable, if I q of the grating cell could be 

reduced , and therefore V qc increased , the efficiency of the grating cell 

would be higher. Table 1 shows that the maximum .1 per unit area obtained 

2 

from a grating cell in our study is 28.8 mA/cm for illumination by a 
tungsten iodide source whose intensity was set to produce the same current 
in the same commercial Si cell as it generated when it was exposed to the AMO 
solar simulator at NASA Lewis Research Center. The tungsten iodide lamp 
is deficient in the blue region when compared with the AMO spectrum so 
that it is possible that I could even be higher than 28.8mA/cm in AMO 

SC 

2 

sunlights. If this same value of I (28.8mA/cm ) had been obtained in the 

* SC 

2 

1-cm grating cell, its maximum power output could be increased by 1.1 times, 

and therefore efficiency would be about 8.8% with AR coating. This indicates 
2 

that the 1-cm grating cell whose characteristics have been measured and shown 
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in Fig. 28 and Fig. 29 is not optimized, higher efficiencies are potentially 

attainable in grating cells. Figure 29 compares the spectral response of 
2 

the 1-cm grating cell to that of the commercial silicon cell whose I-V 
curve is shown in Fig. 28. As expected the grating cell is a blue shifted' 
cell relative to the commerical diffused junction solar cell. 

3.4 Discussion about open circuit voltage 

It is a common practice to represent the I-V characteristic of a solar 
cell under illumination as 

m 

1=2,1 (exp(qV/A kT> - 1) - I (3.1) 

n=l om * * -n .■ • - sc 

where I is a constant with the dimension of current, kT/q is the thermal 
on . 

voltage (0.026 volts at room temperature), the factor A n is a dimensionless 
number having a value larger than 1, I is the light generated current. 

SC 

However, if the dark current flowing in the p-n junction of a silicon 
solar cell consists of two components, viz. 1) the bulk diffusion current, 
and 2) the bulk recombination-generation, current in the transition region, 
then the I-V characteristic of such a cell under illumination can be 
represented approximately in the form 

I = I D0 (-exp(qV/kT) - 1) + i T() (exp.(qV/2kT) - 1)- I gc • (3.2) 

where the first and second terms on the right hand side of this equation 
are bulk diffusion current and bulk transition region recombination current 
respectively, and I DQ , I TQ are constants with the dimensions of current. 

The open circuit voltage V qc can be obtained by solving the equation: 

I D0 (exp(qV oc /kT) - T) + I T() (exp(qV oc /2kT) - 1) = I gc (3.3) 
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for> a conventional silicon solar cell 

If the p-n junctions of a grating type cell have the same characteristics 
as the p-n junction of a corresponding conventional solar cell, then, it is 
obvious , that the of the grating cell is approximately a/b of the 
of the corresponding conventional cell, and I DQ (conventional cell) > I Q0 
(grating cell) > a/b • I q q (conventional cell). Therefore, when the 
I 's of the grating type cell and the corresponding conventional cell 
have the same value, then from Eq. 3.3 and the previous argument, it 
is obvious that the V qc of the grating cell is larger than the V qc of the 
corresponding conventional cell. 

Unfortunately, the experimental results showed otherwise. The ^ oc ' s 

of grating cells made by methods described in section 3.1 are less than 

0.5 volts which is considerably lower than 0.6 volts for convention p/n 

cell and 0.55 volts for conventional n/p cell. Figure 30 shows the dark 

2 

I~V forward characteristics of a 1-cm grating type cell and a commerical 

Si cell with its AR coating removed* It is reasonable to say that the 

excessive large forward current of the grating type cell might be caused 

by some other current component rather than the bulk diffusion current and 

the bulk transition region recombination current* 

(33 34) 

According to Sah 9 , the current flowing in a semiconductor junction 

may be divided into four components according to the location of the recom- 
bination and generation of electrons and holes* These are: 1) the bulk 

recombination -generation or the diffusion current, 2) the bulk recombination 
generation current in the transition region, 3) the surface recombination- 
generation current in the transition region, and 4) the surface channel 
current* The current voltage relationship for these four current components 

may be approximated by I - I exp(qV/mkT) if V > 4kT/q* Analysis shows 

s 

that the reciprocal slope m for the current components 1) to 3) lies between 
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1 and 2, while for the surface channel current component 4) m is greater 
than 2 for silicon junctions and may be larger than 4 for larger channels. 
Figure 31. a is used by Sah to describe the generation recombination locations 
for these four current components. During his investigation, he used a 
field plate to induce a channel on the surface near a junction. A surface 
channel is defined as a region near a p-n junction into which either electrons 
or holes from the other side of the junction can flow without having to climb 
a potential barrier. 

Also shown in Fig. 30, are the reciprocal slopes m for the 1-cm grating 
cell and the commercial cell as a function of forward voltage. For high 
voltages , the large m is due to the series resistance , but around 0 . 2 volts , 
the m for the grating cell is 3.5, and by the interpretation of Sah, it 
means that in this voltage range, the current is dominated by the recombination 
generation current in the channel. Note that for the commercial cell with 
AR coating removed, the m also exceeds 2, this has also been observed by 
Shockley and others 

As to the origin of a surface channel in an alloyed p-n junction grating 
cell, it might be explained as follows. In figure 31. b and figure 31. c 
we sketch the cross-sectional views of a unit composed of a grating line and 
underlaid substrate before and after the alloying process.' The nonuniform 
profile of the p-n junction formed underneath the grating line is a result 
of chemical etching undercut effect on the grating lines. During the alloying 
process , a certain amount of silicon underneath the A1 grating lines would 
mix with A1 (e.g., the eutectic composition of Si and A1 is 11% of Si and 
89% of A1 by weight), and the quantity of Si participating in this process 
is proportional to the thickness of the overlaid A1 film. The p-n junction 
would form at the interface between the molten region and the solid Si 
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substrate. As a result, a very thin layer of heavily doped p+ region 
(about 10 - 100A, the same magnitude as the width of an inversion layer) 
is formed along the edges of the grating line, it is a hole channel. Holes 
from the neighboring bulk p+ region would flow into this channel, a voltage 
drop would develop along the path of current flow due to the high resis- 
tivity of the channel. Part of the hole current would be used to supply 
the recombination current in the surface channel transition region, and might 
be set equal to the change of the current flowing along the channel. The 
recombination current in the surface channel transition region is voltage 
dependent, and is proportional to Sinh v/2, where v is the voltage difference 
across the surface channel transition region. Then, according to the 
analysis by Sah , the channel current may be proportional to Sinh(qV/4kT), 
where V is the voltage difference between the bulk p+ region and the 
substrate . 

I 

Shunting paths also could be responsible for the excessive large current 

(38 ) 

in a grating type cell. Stirn investigated the shunting paths in con- 
ventional cells, and concluded that the shunting paths were localized in 
nature,, the most likely cause of these shunting paths would be scratches 
penetrating through the junction region since in commerical solar cells, 
junction is only 0.3 to 0.5 microns from the surface. For alloyed 
grating type cells studied in this investigation, extremely shallow junctions 
would be formed along the edges of the grating lines for reasons described 
in last paragraph; therefore, shunting paths would be formed most likely 
along the edges of grating lines., 

A possible remedy to eliminate the excessive large forward current in 
a grating type cell would be to use ion implantation technique to make such a 
cell. A1 atoms could be ion implanted into a n-type Si substrate through 
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a suitable mask to define a grating structure p-n junction, and by decreasing 
the energy of the ion beam, an A1 grating structure ohmic contact overlapping 
the p-n junction could be also imposed. In so doing, because there is no 
lateral • spreading of the dopant , no surface channel on extremely shallow 
junction could be formed through metullurgical interaction. 
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4. CONCLUDING REMARKS 

a) The design principles and the underlying theory for the grating-type 
photovoltaic cell have been established. 

b) The grating-type Si photovoltaic cell can have a comparable or even 
greater short-circuit current, and, at least, in theory, a higher open- 
circuit voltage than the conventional Si photovoltaic cell, it is a 
potentially higher efficiency cell. 

c) The grating type cell is a blue-shifted cell relative to the conventional 
cell, it has higher response to blue -ultraviolet light and lower 
response to red-infrared light. 

d) The low open-circuit voltage of the alloyed Al-Si p-n junction grating 
type cells studied in this investigation might be due to the formation 
of extremely shallow junction along the edges of grating lines. 
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5. APPENDIX 

By using the five point formula and the integration method'"'''' of the 
finite difference approximation technique, the second order, self-adjoint 
elliptical partial differential equation: 


-(P(x,y)W x > x - (P(x,y)W y )y + a (x,y)W(x,y) = f(x,y) (5.1) 

can be replaced by a discrete model. The approximate value of W(x^,y^)-, 

viz., w(x.,y.), can be represented in terms of w(x. . y. ), tc(x. . ,y . ) , 

1 1 J ±T± J 

“ (x i sy j-i)’ u(x i»yj + i^» p(x i » y j+i/2** p(x i* y j-i/2 ) ’ 

P(Xi + iy 2 > y j) and P(x._jy 2 ,yj). The positions and some of the symbols are 

shown in Fig. 32. For simplicity, w(x^,y^.) is designated as ok ^ , ‘so are the, 

other symbols. We have 


where 


(i) 


. , — (-L. .to. , »*fR. . .+T.. .0). . -j +B . *01. * . tS • • ) (5.2) 

1,3 D. . 1,3 1-1,3 ■ 1,3 1+1 > 3 1,3 1,3+1 i.3'U,3“l i»3 . 

i,3 


L. . = P. 

i _ 


K. ,+K. 


,j ~ *1-1/2 ,j '2h i _ 1 ’ R i,j, = P i+l/2,.j 2hi * 


K. ,+K. 


T. . = P. 


h . . +h . 

i-l- i. 


h . , +h 
i-l 1 


l i,j ’ *i,j+l/2 2K. ’ B i,j = P i,j-l/2 '2Kj 


h^i+hi k . _ 1 +k. 

D. . = L. . + R. . + T. . + B, . + a. . -±-= — - . 1 . , 

1,3 . 1,3 1,3 1,3 i,3 1,3. .2 ^ 


S. . = f. . 
i,3 i,3 


h . . +h . K . , +K . 

i-l i 3-1 3 


2 . ' 2 

For the successive over-relaxation iterative scheme. 


(m+1) _ . . (m) f (m+1) ^ c (mV 

oil i '- (1— ft) oil • + ft(L. .(j. . . + R. .(j. . 1 

i,3 i,3 1,3 1-1,3 1,3 1+1,3 


+ T. .w. . (m) +l + B. .to. . (m+1) + S. . )/D . '. 

i,3. 1,3 i,3 1,3-1 i,3 i,3 


(5.3) 
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( m+1 ) 

where ok ^ represents the approximate value of VL ^ after m+1 iterations, 

A is the over-relaxation factor, 1.9 was adopted as the value of ft in the 
computation . 

For the computational model as shown Fig., 3, p(x,y) = 1, 
o(x,y) = -1/L n ^ (or -l/l/), f(x,y) = N Q ae /D. 

A mesh is imposed upon the unit cell as shown in Fig. 33 (for the case 
a = 2.pm, b = 10 pm). Note that the mesh spacings are finer near the junction 
and the front light receiving surface, because the variation of carriers 
concentration there is more rapid than in other areas. 

The approximate value of the excess minority carriers at each mesh 
point in the skin region is assigned by U. .(i from 1 to 11, j from 1 to 14); 
and these in the base region is assigned by U. . (i from 11 to 17, j from 1 

3 . 

to 13; and for i from 1 to 17, j from 14 to 21) and V. .. (i from 2 to 41, 

1 » J 

j from 1 to 4 and for i from 1 to 41, j from 5 to 13). 

Because of the boundary conditions on the junction (n'=p'=0) skin region 
and base region constitute two independent systems, therefore the problem 
can be solved separately. 

The general computation scheme is like this ; first , an initial value 
(arbitrary, but a good guess can prompt the convergent rate) is assigned 
to U. . (and V. .), then the iterative procedure according to Eq. 5.3 is 
started. The iterative step is from left to right and from bottom to top. 

After each complete iteration cycle, a test is performed, to see if the 
value U. . (and V. . ) meet the set criterion at each mesh point: 


u. .°’ +1> - 0. . <m) 

_i j2 Lii 

U i,j CllHi) 


(and 


V. < mtl) - V. . 
J-O IjD 


(m) 


v. .(m+1) 
1 


) < e (5.4) 


If the answer is no, then the iteration process is continued, if the 
answer is yes, then it is proceeded to calculate the collection efficiency. 
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The calculation is done according to the - following scheme: 


scs 


skin qNo b 


'‘base 


J scb 
qNo b 


(5.5) 


where I is the short circuit current contributed by the skin region, 

SCS 

I S cb ‘that contributed by the base region. No is the number of photons 
incident onto the cell, b is the half line spacing described in section 2.3. 
And: 


I = -qD 
scs ^ n 


f V d' i dn 
3 junction 


= qo < 


n 


i 10 
K 13 1=1 h i 


U i,13 +U i+l,13 1 

^ ’ ♦ v 


13 

£ 


k. 


U in .+IL . 
10,3 10 , 3+1 


10 j 


=1 ^ 


(5.6a) 


scb 


= q D^ 


f ' . dn 
3 junction 


. 10 U. , C +U. A1 

- ,Dp ( i— I h £ 1-15 2 1<-1 ’ 15 ♦ | . 

14 i=l 1 11 j 3 


U 12,j* U 12,j-H 


Z k. 


(5.6b) 


where n is a unit vector normal to the surface of the region considered. 


After the collection efficiency and the new values of U. . (and V. .) 

i,3 i,3 


are recorded, the whole process is repeated for several times with smaller 
e's, until is is assured that the final result is a reasonably satisfactory 
one . 


3 —1 —3 —3 

For example, for a = 2x10 cm , the e was set at 3x10 , 1x10 , 

—4 —4 

5x10 , and 3x10 in this order and Q base = 1.016 (after 72 iterations), 

0.9521 (after 64 iterations), 0.9165 (after 137 iterations), and 0.9006 
(after 196 iterations) correspondingly, it was deemed satisfactory, and 
the value 0.90 was taken as the approximated value of Qj-, ase * 

The following is a complete computer program used in the computation 


described above. 
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0 IMEhlS ION U< 17,33) , HU (17,33) ,H ( 17) , HK ( 33) 

1HL ( 17, 33) , R ( 17, 33 ) , T( 17,3 3) , D (17, 33) , S (17 ,33 ) , 
1B{I7»33),HL1(41»13),R1( 41 , 13} , Tl( 4-1,13) ,01(41,13) , SI (41, 13) , 
181(41,13) ,V(41 ,13) ,HV(4l, 13) ,G(41),GK(13j ,.f (41 ) ,F1(41 ) , 
1A2(3),A2(3), APHA ( 3) ,N (3 ) 

REAP (5,100) KMAX,ERMAX 

READ (5,600) (A2(L),L=1,3),IA2(L),L=1,3), 
1(APHA(L),L=1,3),(N(L),L=1,3) 

READ (4) ( (U(I;J) ,1=1,17) ,J=1,33) ,( (V(I,J) , 1=1 , 41 ) , J=1 >13) 

L=3 

DO 19 1=1,12 

19 H (I } =0 ♦ 1 
H (13 ) = 0 .2 
HU4) = 0.2 
H (15 ) = 0 .4 
H(16)=0.5 
DO 20 3=1,7 

20 HK(J)=0.2 

DO 21 J=6,15 

21 HK{ J }=0 . 1 
HKt 16)=0. 2 
HK(17) = 0.2 
HK(18)=0.4 

DO 55 J=1 9 , 32 
55 HK( J)=0 .5 
DO 22 1=1, 10 

22 G(I>*1. 

DD 23 1=11,20 

23 G(l)=2. 

DO 24 1=21 , 3fc 

24 G(I)=4. 

G(37>=2. 

G(38) = 2. 

G(39 )=1. 

G140) = 1. 

GKl 1 )=0 ,4 
GK(2 ) =0.6 
GK(3)=0.5 
GK(4)=0 . 5 

DO 25 J=5 , 1 2 

25 GK(J)=1. 

F(l)=APHA(L)*lO.**4l 
DO 31 1=1,15 

31 F(I+ 1 )=F( r)*EXP(-APHA(L)*H(I ) ) 

DD 41 J=2, 13 
DO 41 1=2,1 0 

HL(I,J)=(HK( J-l)+Hk(J))/(2.*H(I-l)) 

R(I, J)=(HK( J-1)+HK(J) >/(2**HCI)) 

T(I, J)=(H(I~1)+H(I))/(2.*HK( J)) 

B(I, J)=(H(Ir 1)+H(I) ) /(2 .*HK( J-l )) 

P(£,J)=HL(I iJ)+R(I,J)tT(I , J )+B ( I, J) 
l+(HlI-l)+HtI))*(HK( J-D+HKIJ) 

42 S (I, J) = (H( 1-1) +H( I) mHK( J-l)+HK(J))'7F(I)/4i 





00 32 J=l,l 3 

32 U(1,J) =10. *0(2,0/11. 

DO 33 1=1,10 

33 U ( 1, 1 ) = l)( 1 , 2 ) 

DO 34 J=2,13 
DO 34 1=2,10 
HU( 1,J J=U( 1,0 

34 U(I, J)=I.9*(S(I,J )+HL(I, 0*U(I-1, J> +R (I , J)*U(I+1, J) 

1,+ T ( I , J )*U CI,J + 1J+BH , 0,*U( I » J-'l) ) /D ( I, 0+ (l »~i • 9} *U'( I, J ) 
DO 35 J=2, 13 
DO 35 1=2,10 

ER=AB5 (U( I, J )-HU(I, O )/U{ X, J) 

IF (ER-ERMAX) 35,39,39 

35 CONTINUE 

Write 16,200 jc,(uit»i3),i=i,io) 

WRITE (6,200) K,(U(l0,O,J=l ,13) 

TU=0. 

00 37 1=1,10 

37 TU=TU+(U(I,l3)+U(I+l,13))*H{I)/2. 

DO 39 J=1 , 1 3 

3ft TU=TU+(U( 10, 0 +U( 10 ,0+1 ) )*HK( 0/2. 
TUP=TU*10./(10.*10.**4l) 

WRITE (6,30 0 K , TUP 
00 TO 50 

39 IF IK .LE. Kf-fAX ) GO TO 30 
WRITE (6,200) KMAX 
50 F(l)=APHA(L)*l0.**4l 
DO 67 1=1,15 

67 F(I + 1)=F( I)tsEXP(-ARHA(U)*H(I )) 

DO 42 J=2, 14 


DO 42 1=12,16 

HL(I,0 = (HMJ-1)+HKU))/(2.*H( I-D) 

R(I, J) = (HK( J~J)+HK(J) )/<2.*H(I) ) 
T(I,J)*(H(I-1)+H(I) )/(2.*HK{ J) ) 

B(I, J)=(H(I-l)+H(I) J/12.*HK(J-1J) 
D(I,J)=HLlI,0+R(I, J)+T(I, 0+5(1, J) 

1 + (H(1-1)+H( D)*(HK( J-1)+HK( 0)/(4‘.*300.*300. ) 

42 S(I,0 = CH(I-1)+H(I) )*(HK< J-1)+HK( J) )*F (l)/4. 
DO 43 J=15, 32 

DO 43 1=2,16 

HUI,J)=(HKU-l)+HklJ))/(2.*H( I-l>) 
R(I,J)=(HK(J-1)+HK(J))/(2.«H(I) ) 

T(I,J) = IH(I-1)+H(I) ) /(2.*HKt J) ) 

B(I,J) = (H(I-1)+H(I) >/(2.*HK< J-D) 

D ( I , J ) =HL (I»0+R(I,J)+T(I,0 + B(I»J) 
1+(H(I-1)+H( I))*(HK< J-D+HKIO) / (4. *300. *300.) 

43 S(I,J)=<H!I-1)+H(I))*(HK( J-1)+HK(J) }*FlI)/4\ • 
F1(1) = APHA(U*10.**41 

M=N(U 
DO 71 1=1, M 

71 Fl( 1+1 )=F1( I )*EXP(-APHA( L)*G(I ) ) 

DO 72 I=M , 39 

72 Fl(I+2) = 0. 
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DO 44 J=2, 12 
DO 44 1=4,40 

HUtI, J)=(GMJ-1)+GKC J)) / (2.*G(I-1) ) 

Rl( I , J )={ GKU-l)tGK(J))/( 2.*G(I)) 

T1(I,J} = (G( I-1)+G(I) )/C2.*GK(J)) 

Bl( I, J)=(G( I-D+GCI) )/(2.*GK(J-l) ) 

OKI, J)=HL1( Ii J)+RKI,j))+Tl (I, J)+B1 (I, J) 

1+(G(I-1) + G( I))*(GK( ) +GK ( J ) ) / (4-. *3 00. *300 . ) 

44 SKI, J) = (G{I-l)+GCm*CGK{ J-D+GM J))*Fl(I)/4. 

K-0 

60 K=K+1 

DO 61 1=12,16 

61 U( 1, 1) ~U( 1,2) 

• DO 62 J= 1 5 , 3 3 

62 0(1, J)=10.*UC2, J)/10.01 
DO 63 1=4,4 0 

63 Vll, 1) =V( X, 2) 

DO 64 1=4,40 

64 V (I j 13)=V (I, 12) 

DO 65 1=1,16 

65 l) (I, 33)=U ( I, 32) 

0(17 , 1 )=( V(4, 1)+U(16, 1) ) /2. 

U(17,3)=( V{ 4,2 )+U <16,3 ) )/2. 

U(17,2)=(U(17,3)+U(17,1)) /2. 

U(17,6) = ( V( 4, 3) +11(16,6) )/2. 

E=UU7,6)-U(17,3) 

U(17,4)=UU7»3)+0.333*E 
U(17,5)=U(17,3)+0.666*E 
U(17 ,9 )=( V(4,4)+U(16»9)}/2. 

E=U(17,9)-U(17,6) 

U(17,7)=U(17,6)+0,4*E 
U(17,8)=U(17,6)+0.8*E 
U(17,14) = (V(4,S)4U( 16,14) )/2. 

. E=U(17 » 14)-U ( 17 ,9 ) 

U117,10)=U( 17,9)+0.2*E 
U (17, ll)=U( 17,9)+0.4*E 
U(17»12)=U(17,9)+0.6*E 
U (17, 13)=U1 17,9) + 0.8*E 
U(17il9)=(V(4,6)+U(16,19))/2. 

E=Ul 17 ,19')-Ul 17 , 14) 

U(17,15)=U(17,14)+0.l*E 
U(17,16)=U( 17, 14) +0 . 2 J .-£ 

U(17,17)=U(17,14)+0.4*E 
U(17ilB)=U(17,14)+0.6*E 
DO 66 J=l,7 

U(17 ,19+2*J) =( V( 4 , 64 J ) +U( 16,19+2*J> ) /2 . 

66 U(17»19+2*J-1,)={U(17,19+2*J) + U( 17,19+2*J-2) )/2. 

DO 68 J=2, 14 

DO 68 1=12,16 
HU( I,J)=U(1,J) 

68 U(t, J)=1.9*(S(I ,J)+HL(I, J) tRU,j)*UI 1+1, J) 

1+T( I , J K-Ul I , J + l ) + B( 1, J)*U ( I , J-l) ) /P(1,J) + (1.-1.9)*U(1,J) 
DO 69 J=15 , 32 
DO 69 1=2,16 
HU(1,J)=U(I,J) 


originao page k 

OP POOR QUALITY 
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69 U(I, J)=1.9*($<I,J> + HUI, J)*U(I-1, J)+R('I,JJ*U(I+i, J) 

1+T(I, J)*U(I,J+1)+B( I, J)*U(I , J-X j)/D(r, J)+(X.-1.9)*U(I, J) 

DO 13 J=1 , 13 
13 V(41,J)=V(40, J) 

V ( 3 i 1 ) =U( 16,1) 

V(3,2)=U( lb, 3) 

V(3,3)=U{ 16,6) 

V(3,4)=U< 16,9) 

V<3, 5)=U(16,14) 

V(3, 6}— U( 16,19) 

DO 70 J=l,7 

70 V{3 »6+J )=U( 16 1 19+2^J) 

DO 73 J=2,12 

DO 73 1=4,40 
HVl I , J )=V (1 , J ) 

73 Vtl, J)=1.9*(S111,J J+HL1 (I , J) *V( 1-1 , J )+Rl ( I, J)*V (1+1 , J) 

■ 1+Tld, J )4V ( X, J + 1J+B1 ( I, J )*V( I , J-l )) /D1 ( I, JJ +11.-1 .9J*V(I, J) 
DO 74 J=2 , 14 
DO 74 1=12,16 

ER=AB$(U( I, >/U( I, J) 

' IP (ER-ERMAX) 74,80,80 

74 CONTINUE 

DO 75 J=15, 32 
DO 75 1=2,16 

E R=A BS ( U l I,J )-HU(I,J))/l)( I , J ) 

IF (ER-ERMAX) 75,80,80 

75 CONTINUE 

DO 76 J=2, 12 
DO 76 1=4,40 

ER=ABS ( V( I, J )— MV (.I»J))/V(-I»J) 

IF (ER-ERMAX) 76,80,80 

76 CONTINUE 
81 7U=0. 

DO 77 1=2,11 

77 TU=TU+(U(I,15)+UII-1,15) )*H(I-l>/2. 

DO 78 J=2 , 14 

78 TU=TU+(U( 12, J)+U(12,J-l))*HK(J-l)/2. 
TUN=TU*i0./(10.*10.**4U 

WRITE (6,400) K,TUN 

WRITE (6,200) K, ( U( 12, J) , J=1 ,33,2) 

WRITE (6,200) K,(U(I, 15) ,1=1,17) 

WRITE (6,200) K,(U(I,25),I=I,17,2),(V(I,9),I=4,40,2) 

WRITE (6,200) K, ( U( I, 19) ,1=1, 17 ,2), ( V ( 1,6) ,1=4, 40,2) 

WRITE (6,200) K, (Ull, 14), 1=12, 17) ,(V(1, 5) ,1=4,40,2) 

WRITE (6,200) K, (U(t, 6), 1=12, 17), (V(I, 3) ,1=4,40,2) 

WRITE (6,200) K-,(V(11, J) , J=1,Z3) 

WRITE (4) l(U(I,J), 1=1, 17), J=l, 33), ((V(I,J), 1=1,41) ,J=1»13) 
GO TO 90 

80 IF (k . LE . UMAX) GO TO 60 
GD TD 81 


90 CONTINUE 
STOP 

100 FORMAT (110, F10. 4) 

200 FORMAT (110/7512. 4/7E12.4/7E12 .4/7612. 4/7E12. 4) 


a?? NAL PA ^B 
OF POOR QUALITY 


300 FORMAT (I10/5H TUP=,E12.4) 



400 FORMAT (IX0/5H TUN=,E12.4) 

600 FORMAT (3E 10 .2/3E1 0 . 2/3E1 0 . 2/31 10 ) 
END 


(ORIGINAL PAGE IS 
OP POOR QUALHY- 



-45- 


6. REFERENCES 

1. Chapin, D. M. , Fuller, C. S., Pearson, G. L.,, J. Appl. .Phys.'- 25, 676 (1954). 

2. Prince, M. , J. Appl. Phys. 26, 534 (1955). 

3. Wolf, M. , IEEE, the 8th Photovoltaic Specialists Conference Record $ 361 
(1970). 

4. Loferski, J. J. , J. Appl. Phys. 27, 777 (1956). 

5. Cummerow, R. , Phys. Rev., 95_ , 16 (1954). 

6. Cummerow, R. , Phys. Rev., 95_, 561 (1954). 

7. Loferski, J. J., Wysocki, J. J., RCA Review, 22, 38 (1961). 

8. Terman, L. M. , Solid State Electronics 2 * 1 (1961). 

9. .Wolf, M. , Solar Energy 5_, 83 (1961).. 

10. Wolf, M., Proc. IRE, 48, 1246 (1960). 

11. Matlaw, S. L. , Ralph, E., L. , J. Appl. Phys., 30 , 541 (1959). 

12. Lindmayer, J. , Allison, J. , COMSAT Tech. Rev., 3, 1 (1973). 

13. Loferski, J. J., Ranganathan, N., Crisman, E. E., Second Semiannual Report 
on NASA Grant 40-002-093, Division of Engineering, Brown University. 

14. Loferski, J. J., Ranganathan, N. , Crisman, E. E., Third Semiannual Report 
on NSAS Grant NGR 40-002-093, Division of Engineering, Brown University. 

15. Grove, A. S., Fitzgerald, D. J., Solid State Electron 9, 783 (1966). 

16. Shockley, W., Read, W. T. Jr., Phys. Rev. 87, 835 (1952). 

17. Grove, A. S., Physics and Technology of Semiconductor Devices- , (John Wiley 

6 Sons, Inc., 1967). 

18. Shockley, W., Bell Syst. Techi J. , 28 , 435 (1949). 

19. Varga, R. S. , Marrix Iterative Analysis ,, (Prentice-Hall, Inc., 1962). 

20. Wackspress, E. L. , Iterative Solution of Elliptic Systems , (Prentice-Hall 
Inc. , 1966) . 



-46- 


21. Wysocki, J. J., Solar Energy £, 104 (1962). 

22. Dash, W. C. , Newman, R. , Phys. Rev. 99 , 1151 (1955). 

23. Philipp, H. R. , Taft, E. A,, Phys. Rev. Letters £, 13 (1962). 

24. Loferski, J. J., Roessler, B., Crisman, E. E., Chen, L. Y., Kaul, R. , 
Walker, J., Ninth Semiannual Report on NSAS Grant NGR 40-002-093, Division 
of Engineering, Brown University. 

25. Van Der Ziel, A., Solid State Physical Electronics , (Prentice Hall, Xnd. , 
1968). 

26. Turner, M. J. , Rhoderick, E. H., Solid State Electron 11 , 291 (1968). 

27. Lepselter, M. P., Sze, S. M. , Bell Syst. Tech. J. 47, 195 (1968). 

28. Polgar, P., Mouyard, A., Shiner, B., IEEE Trans. Elect. Dev. ED-17 , 725 
(1970). 

29. Hansen, M. , Anderko, A., Constitution of Binary Alloys , (McGraw-Hill, 1958). 

30. Saby, J. S., Dunlop, W. C., Jr., Phys. Rev. 90 , 630 (1953). 

31. Biondi, F. J. , Editor, Transistor Technology , Vol. Ill, Chapt. 7, 

(D. Van Nostrand Company, Inc., 1958). 

32. Trumbore, F. A., Bell Syst. Tech. J. 39 , 205 (1960). 

33. Wolf, M. , Ralph, E. L., IEEE Trans. Elect. Dev. ED- 12 , 470 (1965). 

34. Sah, C. T., Proc. IRE, 49_, 1623 (1961). 

35. Sah, C. T., IRE Trans, Elect. Dev. ED-9 , 94 (1962). 

36. Shockley, W., Queisser, H. J, , J. Appl. Phys. 32 , 510 (1961). 

37. Wolf, M., IEEE Eighth Photovoltaic Specialists Conference, August, 1970, 
Conference Record, pp. 360-371. 

Stirn, R. I., IEEE Ninth Photovoltaic Specialists Conference, May, 1972, 
Conference Record, pp. 72-82. 


38 . 



TABLE I 



I sc (mA/cm 2 ) 

WAFER X 
0.l-?0.3&-cm 
<11 1> 

WAFER H 
2-3.5il-cm 
<I00> 

WAFER in 
5 Si - cm 
< II l> 

WAFER EX 
!0£i-cm 

< m> 

WAFER 3X 
8-l2Xl-cm 
< I00> 

GROUP 

NUMBER 

CELL 

NUMBERS 

2o 

(/im) 

2b 

(/tm) 

1 

— if 


31.8 

9.46 

24.7 

23.3 

25.2 

25.2 

2 

2, 8,14 



8.40 

25.8 

23. 9 

27.0 

27.2 

3 

3, 9,15 


95.3 

5.45 

25.6 

wmm 

26.9 

26.1 

4 

4,10,16 


127.6 

3.97 

27.9 

16.7 

27.2 

28.6 • 

5 



190.8 

3.01 • 

26.8 

1 1.8 

. 27.1 

27.5 

6 

6,12,18 

5.1 

254.4 

3.07 

26.7 

6.98 

28.6 

28.8 

















































- 48 - 


TABLE 2 



I.c 

(mA/c m 2 ) 

WAFER 
1 ~ 2X2 -cm 
< 100 > 

GROUP 

NUMBER' 

CELL 

NUMBERS 

2a 

(/tm) 

2b 

(/Am ) 

1 

— 

.10.2 

31.8 

24.2 

2 

2,8,14 

20.4 

63.6 

22.7 

3 

3,9,15 

30.6 

95.4 

21.2 

4 

4,10,16 

40.8 

127.2 

21.9 

5 

5,11,17 

61 .2 

190.8 

19.6 

6 

6,12,18 

81.6 

254.4 

19.3 

7* 

19,22 

20.4 

63.6 

19.7 

8* 

— 

20,23 

40.8 

127.2 

ma 

9 

21,24 



8.4 


* WITH FOUR HORIZONTAL BARS ACROSS THE GRATING LINES 
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FIG. I ONE DIMENSIONAL MODEL FOR A CONVEN- 
TIONAL P/N PHOTOVOLTAIC CELL. 




FIG. 2 TOP VIEW OF A 50 x 50 MIL 2 SMALL 

AREA GRATING CELL WITH 2a = 5.1 fJL m, 
2b = 63.6/i.m. 


not pilmeq' 


bags blank 
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FIG. 3 COMPUTATIONAL MODEL OF A 
DIFFUSED GRATING CELL. 


= S3 p 


(COLLECTION EFFICIENCY) 



FIG, 4 BASE COLLECTION EFFICIENCY AS OBTAINED FROM A NUMERICAL 
AND THE EXACT SOLUTIONS OF A ONE DIMENSIONAL CELL. 



( RELATIVE MAGNITUDE ) 



FIG.5 EXCESS MINORITY CHARGE CARRIER 

DISTRIBUTION IN THE BASE REGION PER 
UNIT VOLUME PER INCIDENT PHOTON OF 
THE ONE DIMENSIONAL CELL AS 
SPECIFIED IN FIG. 4 
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- EXACT SOLUTION 

NUMERICAL SOLUTION 

x Lp = I /u. 




FIG. 6 SKIN COLLECTION EFFICIENCY AS OBTAINED FROM A 

NUMERICAL AND THE EXACT SOLUTIONS OF A ONE DIMEN- 
SIONAL CELL AS SPECIFIED IN FIG. 4. 
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Q (ARBITRARY UNIT) 



FIG. 9 A QUALITATIVE DESCRIPTION OF THE COLLECTION 

EFFICIENCY FOR GRATING CELLS WITH A FIXED VALUE 
OF b/g AND DIFFERENT VALUES OF a 




(^ARBITRARY. UNIT) 



FIG. 10 A QUALITATIVE DESCRIPTION OF THE COLLECTION EFFICIENCY 
FOR GRATING CELLS WITH A FIXED VALUE OF a AND 
DIFFERENT VALUES OF b. 








(X) (cm 2 -sec“ / ti )" 1 


-B0~ 



FIG. 12 NUMBER OF PHOTONS AVAILABLE 

OUTSIDE THE EARTH'S ATMOSPHERE 
AS A FUNCTION OF PHOTON WAVE- 
LENGTH. 



Q (COLLECTION EFFICIENCY) 


B x "GRATING" CELL S 2 =0, S 3 =0(BACK SIDE) 
AvONE DIMENSIONAL CELL S = 0 AT BACK SIDE 



FIG. 13 COLLECTION EFFICIENCY OF GRATING CELL AND A CORRESPONDING 
CONVENTIONAL CELL ' AS A FUNCTION OF PHOTON WAVELENGTH 




io' I0 2 I0 3 . I0 4 I0 5 I0 6 


a(cm“M 

FIG.I4 COLLECTION EFFICIENCY OF A GRATING CELL AND A 

CORRESPONDING CONVENTIONAL CELL AS A FUNCTION 
OF THE OPTICAL ABSORPTION CONSTANT 






Q (ARBITRARY UNIT) 




FIG. 16 A QUALITATIVE DESCRIPTION OF THE BASE COLLECTION 
EFFICIENCY OF GRATING CELLS WITH OPAQUE GRATING 
LINES, b/a = CONSTANT. 
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b 3 »b 2 >b,>b 0 



FIG. 17 A QUALITATIVE DESCRIP 
EFFICIENCY OF GRATING 
LIN.ES. a = CONSTANT. 
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ION OFTHE BASE COLLECTION 
CELLS WITH OPAQUE GRATING 




FIG. 18 THE ALUMINUM -SILICON SYSTEM (FROM 

M. HANSEN, AND A. ANDERKO, CONSTITUTION 
OF BINARY ALLOYS, McGRAW-HILL, NEW 
YORK, 1958.) 
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FIG. 19 TEST SAMPLE CONFIGURATION.' 








FIG. 21 CONFIGURATION OF CELLS PRODUCED ON A 
SINGLE WAFER. 




Q (NORMALIZED QUANTUM 
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FIG. 22 NORMALIZED SPECTRAL RESPONSES ■ FOR CELLS ON WAFER M 
AND A COMMERCIAL SOLAR CELL WITH AR COATING REMOVED 
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0.4 0.6 0.8 1.0 1.2 


WAVELENGTH IN MICRON 

FIG. 23 RELATIVE SPECTRAL RESPONSES 

FOR CELLS ON WAFER I. 2a = 5.lum, 
Isc (roo*) 3 9,46 mA/cm. 



RELATIVE QUANTUM EFFICIENCY 
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WAVELENGTH IN MICRON 

FIG. 24 RELATIVE SPECTRAL RESPONSES FOR 
CELLS ON WAFER m. 2a = 5.l//.m, 

Isc (max) = 23.9 mA/cm 2 . 



RELATIVE QUANTUM EFFICIENCY 
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FIG.25 RELATIVE SPECTRAL RESPONSES 

FOR CELLS ON WAFER X. 2a=5.l/um, 
I sc (max) « 28.8mA/cm 2 ; 



RELATIVE QUANTUM EFFICIENCY 
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WAVELENGTH IN MICRON 

* / 

FIG. 26 RELATIVE SPECTRAL RESPONSES FOR 
CELLS ON WAFER 3ZE. d:b = 8:25 



or 



FIG. 27 A I - cm 2 GRATING 
CELL, 2a = 7. 5^m 


STRUCTURE 
, 2b = 127.2/im. 
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I ( m A /cm 



D .1 . .2 .3 .4 .5 .6 


V (volts) 

FIG. 28 I-V CHARACTERISTIC OF A I- cm 2 GRATING CELL 
UNDER SIMULATED AMD ILLUMINATION 
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CHANNEL NEAR A P-.N JUNCTION. 


I 

'O 

lO 

I 





-8 



3 33 




*5 ^ 

h h y 

> 5 

*w 

7 £/ — 4.0-c 

US 

m 


✓ 

✓ 

11 

✓ 

✓ 

J_L 

_LL 

0 

_la! 

f 

t 

LI 









i 

i 

i 

i 

^rc/ o.o 

6 19 3.0 C 

• / c? pc . 





1 


i 

i 





1 


i 


/a c. .o 



1 

II 

III 

■ 

■ 



1 


II 

II 

III 

■ 

■ 


c. . c- r 

5 2.0 6 

III 


II 

III 

IIJI 

!■ 



III 

■ 1 

■1 

!■! 


mm 





IBI 

■ HI 

■ 1 

■ i 

■■ 

Ill 

hi 

mm m 

■■ 





** 

■ ■■ 
!■! 

■i 

■ 1 

■1 

■1 

111 

■■ 

HI 

III 

SB 

— ”■ 



4 S 1.5 6 

Ill 

III 

■ 1 

■ 1 

■i 

■l 

■■ 
■ ■ 

nft 


— 

— 


III 

II 

III 

II 

III 


■ 


3 6 1.0 <5 

III 

II 

III 

II 

IIJI 


■ 


III 

II 

II 

1 

IfcV 


■ 


, 

III 

II 

II 

1 

ill 




2 3 0.4 0 

III 

11 

II 

1 

llj 




III 

II 

II 

1 

Hi 


■ 


, 2 , 1 

III 

II 

II 

1 

III 


■ 



x (/Am) 0 11.2 1.6 

I(U(I,J)) ' 1 ' 

KVII.J)) I 2 

FIG. 33 DISCRETILIZ ATION 
a * 2 /im , b=IO/xr 




